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! Abstract 



The statistics of the energy and hehcity fluxes in isotropic turbulence are 
studied using high resolution direct numerical simulation. The scaling expo- 
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' nents of the energy flux agree with those of the transverse velocity structure 
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functions through refined similarity hypothesis, consistent with Kraichnan's 
prediction [Q]. The helicity flux is even more intermittent than the energy 
flux and its scaling exponents are closer to those of the passive scalar. Using 



^ ' Waleffe's helical decomposition, we demonstrate that the existence of positive 



X 

' mean helicity flux inhibits the energy transfer in the negative helical modes 



a non-passive effect. 

PACS numbers: 47.27.Ak,47.27.Gs 
The classical theories of fully-developed turbulence [0] were dominated by the concept of 
the energy cascade to small scales. However, kinetic energy is not the only local conserved 
integral of the inviscid equations of motion, the three-dimensional (3D) incompressible Euler 
equations. Since the classical theories were developed, it was discovered that there is a 
second quadratic invariant, the helicity: 
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H{t) = j di^ u(x,t)-a;(x,t). (1) 

Here u is the velocity field and = Vxu the vorticity field. Nonzero mean values of the 
helicity are now known to occur naturally in a wide variety of geophysical fiows, such as 
hurricanes and tornadoes 0. It was proposed in that, if the large-scales of the flow 
are helical (parity non-invariant), then there should be a joint cascade of both energy and 
helicity to small-scales. In that case, the helicity spectrum as well as the energy spectrum 
should satisfy a —5/3 law in the inertial range: H{k) ~ Ch{5 / ^^^^)k~^/^ . Just as for a 
passive scalar, the spectrum of helicity was predicted to be linearly proportional to its mean 
flux 5. This predicted joint cascade was conflrmed recently in a numerical simulation of 
3D Navier-Stokes [§, which seems to support the idea of a relatively minor role for helicity 
in hydrodynamic turbulence 0. 0|1O| give theoretical support to the joint cascade picture. 
However, helicity cannot act exactly as a passive scalar. It is well-known that presence of 
helicity suppresses the nonlinear transfer and acts to inhibit the energy cascade |^,|ll],|12 



In 3D magnetohydrodynamic turbulence, the presence of mean magnetic helicity may affect 



the inertial-range scaling laws [|T3[. Perhaps the most intriguing evidence of a hidden role 
for helicity in hydrodynamic turbulence is from the so-called "shell models" . In a family of 
such models, the GOY shell models, it has been found numerically in |14] that the scaling 



exponents of the energy flux are nearly identical to those for 3D Navier-Stokes precisely 
for the members of the family which have a helicity-like invariant. The statistics of the 



helicity flux itself have also been studied in the GOY models [|15| and in a related class of 



"helical shell models" |16|. However, so far the statistics of the helicity flux have yet to be 
explored in 3D turbulence. It is the purpose of this Letter to study the statistics of energy 
and helicity fluxes in 3D hydrodynamical turbulence by direct numerical simulations, both 
with and without a nonzero mean helicity. 

We have simulated the Navier-Stokes (NS) equation in a 512^ domain. The kinetic energy 



is forced in the flrst two shells [r7|. To add positive mean helicity into the flow, we rotate the 



real and imaginary parts of the velocity vector Fourier amplitude also in the flrst two shells 
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to be always perpendicular to each other with the same handedness |TT]]. The NS equation 
was solved using a pseudo-spectral parallel code with full dealiasing. The time stepping 
was a second-order Adam-Bashforth method with a modified Euler method for the initial 
time step. A statistical stationary state was achieved after five large-eddy turnover times. 
In Fig. 1 we plot the energy and helicity spectra of this final steady-state, in the case with 
mean helicity input. Both spectra have about a decade and a half where a —5/3 power-law 
holds. In Fig. 2 we show for the same simulation the mean spectral fiuxes of energy and 
helicity as a function of wavenumber, normalized by mean energy dissipation e = z/(|Vu|^) 
and mean helicity dissipation 6 = 2z/(Vu : 'Vuj). There is about a decade of inertial range 
where these fiuxes are constant. 

The importance of a local energy flux for studying intermittency in the 3D energy cascade 
was first emphasized by Kraichnan |l|], who used banded Fourier series to define such a 
quantity. The local fiux II^; a(x) measures transfer of energy into small length scales < A 
at a fixed point x in physical space. Kraichnan proposed a refined similarity hypothesis 
(RSH) which relates the scaling exponents Cp of the velocity structure functions, defined as 
(|A£u|^) ~ i'^p, to the scaling exponents Tp of the energy fiux, defined by (|n£;^A|^) ~ A'^^. 
Here |A^u| is the magnitude of the vector velocity increment. Precisely, the RSH relation 
is Cp = 3 + Tp/3- Equivalently, this relation may be stated as the equality (p = (p, with the 
latter defined by 

(|A-n^,Ar/=^)~A^" (2) 

for L ^ A ^ rj, where L is the integral length and r] the dissipation length. To test this 
relation here we shall use instead a smooth filter to differentiate the large-scale and small- 
scale modes, as in our earlier work [|T^. This is the same method used in the large-eddy 
simulation (LES) modeling scheme [|T^. A low-pass filtered velocity u = Ga * u with scales 
< A removed obeys the equation 

dtu + (u- V)u = f - Vp - V-T (3) 
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in the limit of high Reynolds number, where viscous terms can be neglected. Here f , p are 



the filtered forcing and pressure, respectively, and t = u®u — u(S)uis the turbulent stress, 
or spatial momentum transport induced by the eliminated small-scale turbulence. From 
this equation, a balance equation is easily derived for the local density ca = of the 



large-scale energy [T^: 



5teA(x, t) + V-J£;,a(x, t) = Fe,a{^, t) - Ue,a{^, t) (4) 

in which the current J^^a represents space-transport of the large-scale energy, a = f"U 
is the energy input from the force, and 

nE,A(x, t) = - Vu(x, t) : t(x, t) (5) 



is the energy flux out of the large-scales and into the small-scale modes. See also p0|j21 . 

In the same way we can derive a balance equation for the density = u-uJ of the 
large-scale helicity: 

dthAi^, t) + V-J/f,A(x, t) = Fi^,A(x, t) - nj^,A(x, t) (6) 

See [|^. Here 3h,a is a space transport of large-scale helicity and Fh,a is the input from 
the forcing, while 

niT,A(x, t) = -2VlJ(x, t) : t(x, t) (7) 

is the local helicity flux to small scales < A. A set of scaling exponents corresponding 
to the helicity cascade can be defined by 

HA^UhaI"^') - (8) 

for L > A > r]. 

In Fig. 3 we plot the structure functions of the energy and helicity fluxes that appear on 
the lefthand sides of Eqs.(0),(P) for integer values of p from 1 to 9. 



We use the Extended Self-Similarity (ESS) procedure [^] to extract the scaling exponents 
and for 30 < A/r/ < 80. The results are shown in Fig. 4. Together with these we 



plot the scaling exponents of the transverse velocity structure functions. The transverse 
velocity differences are known to be more intermittent than the longitudinal ones p^j25[] 
and thus must dominate in the structure function for which Kraichnan's RSH was proposed, 
which includes all velocity components. As may be seen from Fig. 4, the scaling exponents 
(p and (p are quite close for each p , in agreement with Kraichnan's RSH p|JT8|]. However, 
the scaling exponents for helicity flux are smaller, < (p, indicating that the helicity flux 
is intrinsically more intermittent than the energy flux. This is consistent with the picture 
of the helicity acting similarly as a passive scalar, since it is well-known that the scaling 



exponents of the scalar are smaller than those of the advecting velocity itself [p6l , p7| . As 



shown in the inset of Fig. 4, the scaling exponents of helicity flux are indeed rather close to 



those of passive scalars p8[. It is worth emphasizing that the relation between the scaling 
exponents of energy flux and helicity flux found here is exactly the opposite of that observed 
in the shell models. In [|15| it was shown numerically that the helicity flux in the GOY 



shell model is less intermittent than the energy flux there, while in |T6[ it was shown for 
the "helical G0Y3 model" that energy and helicity fluxes are equally intermittent. Thus, 
despite the fact that the energy flux statistics of the helicity-conserving shell models are 
very similar to those of 3D Navier-Stokes W^, nevertheless the helicity-fiux statistics of 



these shell models are qualitatively different from those of Navier-Stokes. This is discussed 



further in 22 



Fluctuations of the fluxes H£;^a(x, t) and H//^a(x, t) in the joint cascade can also be 
described by the single-point probability density functions (PDF's), which we have calculated 
at various values of A in the inertial range. In Fig. 5 are plotted these PDF's for A/r] = 64. 
fi represents the mean of the flux and a stands for the standard deviation of the flux. The 



PDF of the energy flux agrees with the results reported earlier in There is an obvious 
skewness with a long tail to the right, indicating the forward cascade of energy preferentially 
to smaller scales. To our knowledge, the result on the helicity flux PDF in Fig. 5 is entirely 
new. In contrast to the PDF of the energy flux, it is nearly symmetric. This is because the 
helicity is indefinite in sign and can be both positive and negative, while kinetic energy is 
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always positive. Therefore, the long tail of helicity flux to the left does not indicate backward 
transfer of positive helicity but instead forward transfer of negative helicity. This argument 



can be made precise by means of the helical decomposition of the velocity field [^. One 
sets u = u+ + with 

u±(x,t) = 5:a±(t)h±(k)e^''-^ (9) 

k 

where hs(k) for s = ± are eigenvectors of the curl: 2kxhs(k) = s|k|hs(k). As shown in 
29| , the + modes carry only positive helicity and the — modes only negative helicity. If we 



likewise decompose vorticity as a; = uj^ + a; , we can define 

n|,A(x,t) = -Vu±(x,t):r(x,t), 

n|,A(x, t) = -2VaJ±(x, t) : t(x, t) (10) 

and then n|; ^ represents flux of energy from the large-scale s-modes into the small scales 
and n|^^ represents the corresponding quantity for helicity, with s = ±. See [ p^ . 

In Fig. 6(a) are plotted the 1-point PDF's of for nonzero mean helicity with the 

inset for zero mean helicity. It may be seen that the PDF of the + mode is skewed to 
the right, while that of the — mode is skewed to the left, indicating that both positive 
and negative helicity prefer to cascade forward to smaller scales. We can not see much 
difference between nonzero mean helicity and zero mean helicity. In Fig. 6(b) are plotted 
the 1-point PDF's of n| ^ for nonzero mean helicity with the inset for zero mean helicity. 
Although the PDF of 11^ ^ is similar to that for the total energy flux itself, and skewed to 
the right, the PDF of 11^ ^ is much more nearly symmetrical for nonzero mean helicity case . 
Therefore, the positive mean flux of helicity seems to inhibit forward transfer of energy in the 
— modes, exhibiting the "non-passive" character of helicity in fully-developed turbulence. 
This discovery is in agreement with the assertion of Kraichnan [|^ that interactions will be 



strongest between the helical modes of the opposite signs. See also p9 |. 

To see further how the helicity flux affects the energy flux we have calculated the PDF of 
the energy flux conditioned on the value of the helicity flux, in the case with 6^0. In Fig. 7 
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are shown two PDF's of energy flux, one conditioned on helicity flux in the range [0.0, 3a^] 
and the other in the range [3a^ , 6a^]. We can see that energy flux is more intermittent 
when the helicity flux is small than when the helicity flux is large. Clearly, the tail of 
the PDF with the large helicity flux is reduced compared to that with the small hehcity 
flux, implying that large helicity flux "blocks" large energy flux. However, the effect of the 
conditioning upon helicity flux is not entirely simple. While very large energy flux events 
are suppressed in the presence of large helicity flux, the probability of larger fluctuations 
near the center of the distribution is actually increased. It is an interesting question whether 
similar effects would be seen on energy flux when conditioning on values of a passive scalar 
flux. The effects observed imply only statistical correlations between two fluxes and no 
direct dynamical interaction. A detailed theory of the phenomena remains to be developed. 
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Figure Captions 

FIG. 1. Energy and helicity spectra. 

FIG. 2. Normalized energy and helicity fluxes versus wavenumber. 

FIG. 3. Structure functions (a) of energy flux, (b) of helicity flux. The scaling exponents 
are measured using the data between two dashed lines. 

FIG. 4. Scaling exponents of transverse velocity increments, energy flux, and helicity flux. 
In the inset we show the scahng exponents of hehcity flux(«) and passive scalar(O) 

FIG. 5. 1-point PDF's of energy and helicity fluxes. 

FIG. 6. 1-point PDF's of (a) ± hehcity fluxes, (b) ± energy fluxes for different mean 
helicity inputs 6 when A/77 = 64. The inset is for 5 = 0.0. 

FIG. 7. Normalized energy flux PDF's conditioned on helicity flux when A/77 = 64, 5 7^ 0. 
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